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We study 2D quantum gravity on spherical topologies using the Regge calculus approach. Our goal is to 
shed new light upon the validity of the Regge approach to quantum gravity, which has recently been questioned 
in the literature. We incorporate an R 2 interaction term and investigate its effect on the value of the string 
susceptibility exponent 7 3 tr using two different finite-size scaling Ansatze. Our results suggest severe shortcomings 
of the methods used so far to determine 7 3 tr and show a possible cure of the problems. To have better control 
over the influence of irregular vertices, we choose besides the almost regular triangulation of the sphere as the 
surface of a cube a random triangulation according to the Voronoi-Delaunay prescription. 



1. INTRODUCTION 

2D Euclidean quantum gravity is believed to 
be an important toy model on our way to a real- 
istic 4D Minkowskian quantum theory of gravity. 
Analytic calculations using conformal field theory 
and matrix models have led to a remarkably good 
understanding of the 2D theory. In higher dimen- 
sions, however, numerical methods, like Regge 
calculus [1], will probably be indispensable tools 
to probe quantum gravity. One would therefore 
like to recover all known continuum results. One 
interesting aspect is the string susceptibility ex- 
ponent 7 s t r , which is the sub-dominant correction 
to the large area behavior of the partition func- 
tion Z(A) oc A~ lstI ~ 3 e~ XRA , where Xr denotes the 
renormalized cosmological constant. The expo- 
nent 7 s t r depends on the genus g of the surface 
through the KPZ formula 7str = 2 - |(1 - g) [2]. 

In Regge calculus one introduces a curvature 
square term and deduces from its expectation 
value an estimate on 7 str through a finite-size scal- 
ing (FSS) analysis. The model for R 2 gravity is 
defined by 

with the gravitational action taken as Sq = 
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J d 2 Xyfg(\ + § R + f R 2 ). The coupling constant 
a sets a length scale of ^fa, and A := A/a can be 
used to distinguish between the cases of weak R 2 - 
gravity (A^> 1), where the KPZ scaling is recov- 
ered, and strong i? 2 -gravity (A <C 1) where it was 
found [3], that Z(A) oc A<^- 3 e~ s -l A e~ XRA - bA , 
with the classical action S c = 167r 2 (l — g) 2 , some 
constant b, and 7g tr = 2 — 2(1 — g). Note that 
only for the torus (g = 1) weak and strong 
R 2 gravity have the same scaling behavior with 

7str = 7str = 2- 

For the sphere (g = 0) and the dl/l measure, 
Gross and Hamber [4] found weak numerical ev- 
idence for 7 s t r = —1/2, while Bock and Vink [5] 
reported a failure of the KPZ formula. However, 
in Ref. [6] it was claimed that the string suscepti- 
bility exponent 7g tr for strong R 2 gravity was con- 
sistent with the theoretical predictions for g = 0, 
but not for g = 2. Due to those incompatible 
results we started to reinvestigate the case of R 2 
gravity on the sphere in detail. 

2. METHOD AND FSS ANSATZE 

We simulated the partition function 

Z = J V^l) exp (- J2 ( XA * + ( 2 ) 

i 

where T>ij,(1) is the most commonly used "com- 
puter" measure dl/l. The notation is identical to 
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that used in Ref . [7] . The deficit angle is given by 
Si = 2ir — ^2 tDi @i(t) with 6i(t) being the dihedral 
angle at vertex i, and Ai are barycentric areas. 
The only dynamical term is the i? 2 -interaction, 
because we held A fixed during the update. As 
global lattice topology we used the triangulated 
surface of a three-dimensional cube [5] , and a ran- 
domly triangulated sphere constructed according 
to the Voronoi-Delaunay procedure. In this way 
we can control the influence of irregular trian- 
gulations. For spherical topologies we have the 
relations N - 2 = N 2 /2, N - 2 = #i/3, and 
2Ni = 3N 2 , where #0,^1, and N 2 denote the 
number of sites, links and triangles, respectively. 

2.1. Finite-size scaling 

The methods of Refs. [4,5] have some intrinsic 
inconsistencies. It was shown in Ref. [5] that the 

dimensionless expectation value R 2 := a(^ 8 - -j-) 
depends only on N 2 and the dimensionless param- 
eter A. Sending N 2 — > 00 one expects R 2 (A,oo) 
to be expandable in a power series whose first 
three terms read as 

R 2 (A) = b A + b 1 + b 2 /A, (3) 

where b = — (Xr - X)a, &i = 7 str - 2 for A > 1, 
and b = -b - (X R - X)a, &i = ~f' stT - 2, b 2 = S c 
for A « 1. 

The scaling Ansatz of Ref. [5] is to consider first 
a power series expansion of R(A, N 2 ) in N 2 . The 
expansion of R 2 is, however, not done at fixed 
A, but at a fixed discretization scale set by the 
dimensionless parameter do := A/N 2 : 

R 2 (a , N 2 ) = N 2 c (d ) + c 1 (d ) + c 2 (a )/N 2 • • -(4) 

A similar good guess would have been to consider 
an expansion in a linear length parameter \/N 2 . 
In a second step the coefficients c 8 are expanded in 

(0) , - (1) , (0) , - (1) , 

a as c = c ' + a c ' + . . ., ci = c\ '+a c\ ' + . . ., 

and c 2 = [c^ + a^c^ + ■ ■ ]/ao, and then the 

continuum limit ao — > is taken. A comparison 

with (3) then yields Cg 1 ' = &o,c^°' 1 = bi,d 2 °^ = 

b 2 . But because ao is fixed, and not A, there is 

no clear control over the crossover from weak R 2 

scaling behavior (A 3> 1) to strong R 2 scaling 

behavior (A 1). 

We therefore suggest an alternative approach, 

where we look at the FSS behavior at a constant 



value of A. Expanding R(A, N 2 ) at constant A 
we obtain 

R{A,N 2 ) = N 2 do{A) + d 1 {A) + d 2 {A)/N 2 + ...{5) 

Also here other FSS Ansatze are possible. The 
next step is to expand the coefficients rf 8 - as a 
power series in A. The coefficient d\ carries all 
the necessary information to extract the string 
susceptibilities. A comparison with (3) yields 
di(A) = b A + 7str - 2 + 0(1/ A) for A > 1 and 
di(i) = S c /A + 1 ' stT -2 + b A + O(A 2 ) for A « 1. 
If we plot di versus A we expect to see a linear 
behavior for very large A, and a divergent behav- 
ior for small A, from which we can extract 7 str as 
well as 7^ tr . 

2.2. Simulation parameters 

To update the links we used a standard multi- 
hit Metropolis update with a hit rate ranging 
from 1 ... 3. We ran simulations at constant val- 
ues of 1/do ranging from 0.2 - 1 280, going very 
far beyond the range of Ref. [5] , where only very 
small values of 1/do in the range of 0.5 - 5 were 
studied. The simulations at constant A were per- 
formed in the range of 14 - 7 800. Usually, the size 
of the lattices varied from 218 up to 17498 lattice 
sites, corresponding to 648 up to 52 488 link de- 
grees of freedom. For each run we recorded about 
10 000 measurements of the curvature square R 2 
on every second MC sweep. The statistical errors 
were computed using standard jack-knife errors 
on the basis of 20 blocks. The integrated autocor- 
relation time Tf, 2 of R 2 was usually in the range 
of 5 - 10. 

3. RESULTS 

3.1. Scaling at fixed do 

For the 4 largest values of do the scaling Ansatz 
(4) seems to works well even without the c 2 coef- 
ficient. A closer look on the curves for 1/do = 10 
(1/do = 20), see Fig. 1, shows apparently two 
scaling regions, divided approximately by a line 
through 1/N 2 ^ 0.0005 (0.0003). We interpret 
this region as the crossover region from A ^> 1 to 
A < 1. Because A was chosen to be N 2 /2a, we 
decrease A either by decreasing N 2 or by increas- 
ing a. This means we always start out on small 
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Table 1 

Fit results for Ansatz Eq. (4) (do = l/2a) 
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Figure 1. Scaling of R 2 /N 2 versus l/N'j for 
l/a = 10 and l/a = 20. 
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Figure 2. c\ versus 1/ao for different FSS 
Ansatze. 



lattices in the strong R 2 gravity regime, and end 
up on sufficiently large lattices always in the weak 
R 2 gravity regime. It is therefore hard to imag- 
ine that one can fit the whole range of data points 
with the same Ansatz (4) . For the lower values of 
fio, the data points of the smaller lattices begin 
to show a clear deviation from the straight line 
behavior, which means, that the c 2 coefficient, 
originating from the exponential damping factor 
of the classical action S c for the case 1, wins 

more and more. To obtain more data points for 
small A we ran these simulations on lattices as 
small as No = 26. We then fitted all data points 
using Ansatz (4). In Refs. [5,6] it is claimed that 
in the limit of do — > 0, c\ should approach 7 str — 2, 
and the value of c 2 do should then approach 167T 2 . 
This is, of course, only partially true, because 
technically this limit is approached by increas- 
ing a, and therefore decreasing A, so that strictly 
speaking c\ should approach 7g tr — 2. Discard- 
ing the measurements on the larger lattices leads 
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to an improvement, because they belong to the 
regime where A is large and the c 2 coefficient is 
not so important. One needs very small lattices 
to expect good results. Then, however, one deals 
with such small lattices, that one needs to worry 
about finite-size effects. 

The values of the fit parameters along with the 
number of degrees of freedom (dof) for an ac- 
ceptable quality of the fit can be found in Ta- 
ble 1. As a final estimate of this analysis we take 
the average of the four MC estimates. We feel 



that y st 



-1.9(3) is a conservative estimate 



that is nevertheless still consistent with the pre- 
diction 7g tr = 0. Another way to improve the 
quality of the fits is to include the next order 
correction term to (4). We added two possible 

1/2 

correction terms, one of the form c 3 / 2 (do)/N 2 
and one of the form 03(0,0) /N$. In Fig. 2 one can 
see the values of the coefficient c\ obtained ac- 
cording to three different FSS Ansatze. For small 
values of do all three curves seem to come closer 
together, and are approximately around the theo- 
retical value of 7g tr . For large values of do one sees 
no way how to extract 7 str . If one extrapolates 
to large values of do, then this means A> 1, but 
in the same limit the discretization scale becomes 
very large. 

3.2. Scaling at fixed A 

Looking at the raw data in the plot R 2 /N 2 ver- 
sus I/N2, see Fig. 3, we first note that all curves 
for different values of A are significantly curved. 
Straight line behavior is visible only asymptoti- 
cally for large values of N 2 ■ 

To determine the coefficient di we used Eq. (5), 
and also a three-parameter fit of the form 

R(A, N 2 ) = N 2 d (A) + d^A) + d 3/2 (A)/N^,(6) 

which yielded a better, but still fairly large total 
\ 2 ■ In Fig. 4 we observe that qualitatively both 
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curves fulfill the theoretical expectations, namely 
they show a divergence at small A and a flatten- 
ing slope at large A. For large A we note that 
no clear linear slope can be observed, which sug- 
gests the presence of FSS corrections to scaling. 
Unfortunately, there are not enough data points 
to obtain a precise estimate for either 7 str or 7g tr . 



R7N 

2 

0.440 





A=1825 


□ 


A=912 





A=456 


X 


A=228 


+ 


A=1 14 


A 


A=57 




A=28 


o 


A=14 



Figure 3. R 2 /N'j versus I/N2 with A = const. 
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Figure 4. di versus A for the randomly triangu- 
lated sphere for two different FSS Ansatze. 



4. CONCLUSION 

We tried to measure the string susceptibility 
exponent 7 str for spherical topologies using two 
FSS Ansatze in R 2 . Using the irregular trian- 
gulated surface of a cube and a randomly trian- 
gulated sphere, we observed no qualitative differ- 
ence between the two types of lattices. The re- 
sults we obtained for the dl/l measure show that 
the Ansatz (4) is not applicable to determine 7 str . 
The extrapolation of do — > can only predict 7g tr , 
but with the caveat in mind that already in the 



first extrapolation (4) one mixes data with small 
and large A. Moreover, to reach a sufficiently low 
value of A, one has to extrapolate to very small 
system sizes. In this way one will necessarily en- 
counter large finite-size corrections. As indicated 
by our result of ~f' stT = 0.1(3), the corrections seem 
to be unimportant for the sphere. However, we 
suggest that they are the source of the "failure" 
for the bi-torus in Ref. [6]. 

In principle the Ansatz to work at a well con- 
trolled A, as in (5), should be capable to predict 
both 7 s t r and 7g tr , but we have experienced large 
FSS corrections to scaling, which make it very 
difficult to extract the coefficient di with high ac- 
curacy. The results we obtained can still be con- 
sistent with the theoretical prediction of Ref. [3] 
for both, 7 s t r and 7g tr . 

In light of these difficulties it seems to us very 
questionable if the method of using FSS proper- 
ties of the expectation value of the R 2 term will 
ever reach a high enough accuracy to probe the 
formulas for 7 str and f' stT . One needs to develop 
more direct approaches to measure 7 str , as has 
been done for the DTRS method [8]. The main 
conclusion is that for the dl/l measure a failure of 
Regge calculus has not yet been shown. Details 
of this study will be presented elsewhere [9]. 
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